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Abstract 

■PT-symmetric Hamiltonians and transfer matrices arise naturally in statistical mechanics. These 
classical and quantum models often require the use of complex or negative weights and thus fall 
outside of the conventional equilibrium statistical mechanics of Hermitian systems. T^T-symmetric 
models form a natural class where the partition function is necessarily real, but not necessarily 
positive. The correlation functions of these models display a much richer set of behaviors than 
Hermitian systems, displaying sinusoidally-modulated exponential decay, as in a dense fluid, or 
even sinusoidal modulation without decay. Classical spin models with VT symmetry include Z{N) 
models with a complex magnetic field, the chiral Potts model and the anisotropic next-nearest- 
neighbor Ising (ANNNl) model. Quantum many-body problems with a non-zero chemical potential 
have a natural PT-symmetric representation related to the sign problem. Two-dimensional QCD 
with heavy quarks at non-zero chemical potential can be solved by diagonalizing an appropriate 
PT-symmetric Hamiltonian. 



I. INTRODUCTION 



The fundamental importance of VT symmetry was first pointed out by Bender and 
Boettcher in their seminal work on quantum- mechanical models [B]. Their work grew out of 
the observation that the Hamiltonian 

H=p^ + igx^. (1) 

has only real eigenvalues. Bender and Boettcher observed that the Hamiltonian H, while 
not Hermitian, is invariant under the simultaneous application of the symmetry operations 
parity V : x ^ —x and time reversal T : i ^ —i. This symmetry ensures that all 
eigenvalues of H are either real or part of a complex pair. The argument is simple: if 
Hlip) = E\^) then HVT |^) =VTH\i>) = VT E\i>) = E*VT \ip). Thus if E is an 
eigenvalue, E* is an eigenvalue as well. The unexpected feature of the ix^ model discovered 
by Bender and Boettcher is that all of its eigenvalue are real. Since that discovery, research 
on PT-symmetric systems has grown enormously, with applications in many areas of physics 

mm- 

The ix^ model has its origins in the study of the Lee- Yang theory of phase transtions 
\'27\ Lee and Yang showed how the critical properties of the Ising model could be 

studied using an analytic continuation of the external magnetic field h to imaginary values. 
They demonstrated that the zeros of the partition function were all on the line Re (h) = 0, 
reaching the real axis at the critical point. It was subsequently realized that the phase 
transition the edge of the gap in the distribution of zeros above the critical temperature is 
itself a critical point described by an field theory [20j. In two dimensions, this field theory 
at its critical point yields the simplest of the non-unitary minimal conformal field theories 
[T2| [T3| HH 05]. In one dimenson, the field theory reduces to ix^ quantum mechanics. 
There is a connection between PT-symmetric quantum mechanical models and conformal 
field theories in two dimensions [19]. This connection is closely related to the proof that 
a large class of PT-symmetric quantum mechanical models, including the ix^ model, have 
only real spectra [T7]; see also [18j. 

Our focus here will be the application of generalized VT symmetry to some well-known 
models of statistical mechanics. By generalized VT symmetry, we mean that the relevant 
unitary operator is not necessarily the parity operator; in many cases, the role of V will be 
played by the charge conjugation operator. Some models, such as the Z{N) spin model with 
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a complex magnetic field, will have a direct connection to Lee- Yang theory. In other cases, 
such as the anisotropic next-nearest-neighbor Ising model (ANNNI model) |39| or the chiral 
Potts model [21t [55] . no complex numbers will appear in the formulation of the model and 
the appearance of VT symmetry will be hidden. 

For those already familiar with the subject of VT symmetry, we should emphasize that 
the focus here is somewhat different from the majority of work on the subject. We are 
interested in extant models in statistical mechanics in which PT symmetry plays a role, and 
not on the statistical mechanics of "PT-symmetric models per se. There are some differences 
between the application of PT symmetry to statistical mechanics and work on PT-symmetric 
quantum mechanics. For the models we discuss, the question of integration contours in the 
complex plane does not arise. Furthermore, the standard inner product, as opposed to the 
CPT inner product often used in PT-symmetric quantum mechanics, is sufficient. Thus 
the models we discuss are similar to the Ising model in an imaginary field, which is the 
prototypical example of VT symmetry in statistical mechanics. 

Our original interest in this subject was motivated by the sign problem, particularly its 
appearance in QCD at finite density |29t HO) . Essentially, the sign problem is really an 
instance of the general problem of complex weights in a statistical sum. Such weights may 
arise because the Hamiltonian or action is complex; as we explain below, the problem arises 
naturally in Euclidean quantum field theories with non-zero chemical potential. The sign 
problem is a large barrier to first-principles lattice simulations of QCD at finite density and 
the study of color superconductivity |T]. As we show in section |V| all Euclidean quantum 
field theories with a sign problem due to a non-zero chemical potential have a generalized 
VT symmetry. The class of statistical models with generalized VT symmetry is very large. 
We will show below that this class is precisely the set of models for which the complex 
weight problem can be reduced to a sign problem, i.e., there is a representation in which 
the partition function is constructed from only real, but not necessarily positive, weights. 
We believe that generalized VT symmetry is likely to play a role in any future solution of 
the sign problem. 

The spectral properties of statistical models with VT symmetry are different from those of 
Hermitian theories, which have the property of spectral positivity. Essentially, this property 
ensures that connected two-point correlation functions of observable will fall monotonically 
to zero with separation. When there is a mass gap between the ground state and the lowest 
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excited state, this mass gap determines a minimum rate of exponential decay for correlation 
functions. Systems with generalized VT symmetry need not have spectral positivity, and 
their two-point correlation functions can show sinusoidally-modulated exponential decay, as 
in a dense fluid, or even sinusoidal modulation without decay, as in periodic phases. Because 
these behaviors are observed in nature, it is perhaps unsurprising that PT symmetry occurs 
in many problems of statistical mechanics. 

The remainder of this article is structured as follows: section II introduces the funda- 
mentals of generahzed VT symmetry and its role in the sign problem. Section III discusses 
the general spectral properties of T'T-symmetric models. In section IV, we treat several 
models from classical statistical mechanics. Section V consider quantum statistical models. 
Wc show how VT symmetry makes it possible to solve a form of two-dimensional QCD at 
finite density. A final section provides brief concluding remarks. 

II. FUNDAMENTALS OF VT SYMMETRY AND THE SIGN PROBLEM 

In this section, we discuss some fundamental aspects of VT symmetry. We also discuss 
here some fundamental aspects of the sign problem that have a close relation with VT 
symmetry. Although we will for the most part consider PT-symmetric Hamiltonians, in 
classical statistical mechanics it is often most natural to work with PT-symmetric transfer 
matrices, and the results for Hamiltonians have obvious counterparts for transfer matrices. 
Typically, the transfer matrix T of a lattice model acts as a Euclidean propagator similar to 
exp {—tH) in a particular lattice direction. If the lattice has length L in that direction, then 
with periodic boundary conditions the partition function is Z = Tr (T^) . It is sometimes 
convenient to identify L with (5 and T with exp {-H) so that Tr (T^) can be written as 



A. Eigenvalues of T^T-symmetric systems 

Given a Hilbert space, the adjoint of an operator H is defined using the inner product: 




(2) 



which we often write as 



*T 



(3) 
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Operators satisfying = H are said to be Hermitian (we follow the usual physics practice 
and do not distinguish between Hermitian and self-adjoint operators). It is a standard 
result that the eigenvalues of a Hermitian operator are real. Operators with a generalized 
VT symmetry have a more general constraint on their eigenvalues: they are either real 
or form a complex conjugate pair. A Hamiltonian H is PT-symmetric if [H^VT] = 0. 
Let lip) be an eigenstate of H with eigenvalue E. Then we have HVT =VT H = 
VT E = E*VT Thus we see that VT is an eigenstate of H with eigenvalue 
E*. If is an eigenstate of VT , then necessarily E is real. However, VT- symmetric 
Hamiltonians can also have complex conjugate pairs of eigenvalues. The case where two or 
more eigenvalues are not real is usually described as broken VT symmetry. Eigenstates of 
H associated with a complex eigenvalue pair cannot be eigenstates of VT. In the models 
we consider, T is implemented as complex conjugation, and P is a unitary operator obeying 
[V,T]=0 and p2 = 1. 

B. Bender-Mannheim theorem 

Not all PT-symmetric models are obviously so. A simple criterion for a PT-symmetric 
Hamiltonian H (or transfer matrix T) has been given by Bender and Mannheim pjOJ. If 
the characteristic polynomial det [H — XI] has real coefficients, then H has a generalized 
VT symmetry. Interesting models arise in statistical mechanics with hidden VT symmetry 
when the transfer matrix T is real but not symmetric, as will be discussed below in section 

m 

The striking feature of the ix^ models is that all of its energy eigenvalues are real. This is 
usually referred to as unbroken VT symmetry. We will distinguish three different behaviors 
of "PT-symmetric models, based on the reality of the eigenvalues of the Hamiltonian H or 
transfer matrix T. We order the eigenvalues of H by their real parts, so that the ground state 
of H has the eigenvalue with the lowest real part. Similarly, we order the eigenvalues of T by 
their magnitude; the eigenvalue largest in magnitude is the analog of the ground state energy. 
Typically the different behaviors are each associated with a different part of parameter space. 
In region I, VT symmetry is unbroken, and all eigenvalues are real. In region II, the lowest 
eigenvalue of H is real, but VT symmetry is broken by one or more pairs of excited states 
becoming complex. For a transfer matrix T, the eigenvalue largest in absolute value is real 
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in region II, but other eigenvalues are complex. In region III, VT symmetry is broken by 
the ground state of H becoming complex. For a transfer matrix T, the eigenvalue of largest 
absolute value becomes complex. Strictly speaking, by our definition above this means there 
are two ground states. This does not occur in conventional quantum-mechanical systems 
with a finite number of degrees of freedom, but does occur in PT-symmetric models. It also 
occurs, of course, in Hermitian models with an infinite number of degrees of freedom, and is 
the basis of spontaneous symmetry breaking. As we show in the next section, this distinction 
is physical, and manifests directly in correlation functions. In region II, oscillatory behaviors 
appears in correlation functions. In region III, the system is in a spatially modulated phase. 
We emphasize that the behavior of correlation functions seen in regions II and III cannot be 
obtained from conventional models for which H is Hermitian: such behavior is incompatible 
with the spectral representation of the correlation function for Hermitian theories. 



C. Connection to complex weight problem 

A naive treatment of PT-symmetric models often involves a sum over complex weights, 
as in the path integral treatment of the ix^ model or in the continuation of the Ising model 
to imaginary magnetic field. In many areas of physics, there are problems where we wish to 
evaluate the "average" of some quantity x over an ensemble with complex weights Wj eV- : 

{x) = 4 

There is no effective general algorithm for calculating such sums, as is the case for positive 
weights; see, e.g., [HI [T5| |28]. This problem is generally referred to as the sign problem, be- 
cause even the case of negative weights is difficult. We will discuss the sign problem in detail 



for several models in sections IV and |V| but here discuss the general role of PTsymmetry. 
Let us consider the case where the w/s can be written as exp {—PEj), as is typical in sta- 
tistical mechanics where Ej is eigenvalue of some operator H and /3 is the inverse of the 
temperature T. For a Hermitian system, the partition function 

Zm = J2e-^^^ (5) 
j 

is real and positive for all real values of /3, because the eigenvalues are all real. On the other 
hand, if H is PT-symmetric, the eigenvalues are either real or occur in complex conjugate 
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pairs. We will prove in sections III that Z may be represented as 



Z = 5^ e-^^^ + {e-^^^ + e-^^.*) (6) 

real pairs 

which is real but not necessarily positive. As a consequence of the Bender-Mannheim the- 
orem, we have the following characterization of partition functions: Models which are Her- 
mitian, or equivalent to Hermitian models under a similarity transform, have Z real and 
positive for all real /3; models which have generalized VT symmetry have Z real, but not 
necessarily positive, for all real /3. A Laplace transform argument shows that the converse 
is also true. Thus models with generalized VT symmetry are precisely the class of models 
in which the complex weight problem can be reduced a genuine sign problem, the problem 
of averaging over positive and negative weights. 

D. Equivalence to Hermitian if VT symmetry is unbroken 



The difficulty presented by the sign problem depends directly on VT symmetry breaking 
or its absence. Mostafazadeh pO] has proven that when VT symmetry is unbroken (region 
I) and the spectrum is non-degenerate, there is a similarity transformation 5* that trans- 
forms a PT-symmetric Hamiltonian H into an isospectral Hermitian Hamiltonian via 
Hh = SHS^^. If 5* can be found, the transformation of H into Hh eliminates the sign 
problem for PT-symmetric quantum Hamiltonians throughout the interior of region I. On 
the boundary of region I, where two or more eigenvalues become degenerate, it is possible 
that the Hamiltonian can be of a non-diagonalizable Jordan block form [9j. The equiva- 
lence to a Hermitian Hamiltonian also applies to PT-symmetric transfer matrices T, but a 
further restriction to positive eigenvalues for T is necessary for the elimination of the sign 
problem. Thus there are regions of parameter space where the sign problem can be removed 
by a similarity transformation. Unfortunately, the explicit construction of the similarity 
transform typically requires knowledge of the exact eigenvalues and eigenvectors. However, 
there are some models, notably the — Ax^ model, for which the similarity transform or an 
equivalent functional integral transformation is known |2[ E [251 ES]. In regions II and HI, 
the sign problem has an underlying physical basis, and cannot be removed by a similarity 
transformation. The negative weight contributions to the partition function Z arise from the 
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contributions of complex conjugate eigenvalue pairs associated with VT symmetry breaking. 
It is that breaking that in turn gives rise to the oscillatory and damped oscillatory behavior 
of two-point functions characteristic of many physical systems. 

E. Real Representations of PT-symmetric Hamiltonians 

There are VT symmetric models like the ANNNI model |39| where the classical Hamilto- 
nian, corresponding to the action in the path integral formalism, is real. Such a model can 
be simulated with no difficulties of principle throughout its parameter space. In Hermitian 
systems, the existence of an antiunitary involution commuting with the Hamiltonian implies 
that there is a basis in which H is real; see, e.g., [22j. This theorem easily extends to the 
case of those PT-symmetric systems for which (PT)^ = 1, and can be applied to transfer 
matrices as well as Hamiltonians. This suggests the existence of a class of PT-symmetric 
models which can be simulated in all three regions, but no general criterion for determining 
the class is known. 

The proof that for any PT-symmetric Hamiltonian there is a basis in which the matrix 
elements of H are real follows closely the proof for Hermitian systems |22]. The antiunitary 
operator VT commutes with the Hamiltonian \PT,H] = and satisfies (VT)^ = 1. These 
properties are sufficient to construct a real representation of the Hamiltonian H. The first 
step is the construction of a PT-invariant basis ipa- We start from any non-zero vector (pi 
and complex number ai. The vector ipi = ai0i +PTai0i is invariant under VT. Choosing 
a vector 02 orthogonal to ipi, we form the vector tp2 = ^202 + VTa2(f)2, where a2- The inner 
product of ip2 and ipi is zero: 

(Hi^i) = a*2 (02|^i) + "2 {VT<j)2\^Pi) = a*2 {{VTf 02|^i> + «2 (Pr^2|V'i) (7) 

(V'2|V'l) = 01*2 (02|V^l) + «2 (Pr02|^l) 

= o + «2((pr)'02|pr^i>* 

= "2 (02|V'l)* = 

The complex number a2 can be adjusted to ensure that ip2 is not the zero vector. By 
proceeding in this fashion, a PT-symmetric orthogonal basis can be constructed. In this 



8 



basis, the matrix elements of the Hamiltonian are reah 



= {VTtlJalVTH^,)* 
= {VT^a\H{Vr)Ar 

= {^aWb)* 

TT* 

= ^ab 

Note that this results holds independent of the appearance of complex eigenvalues in H. 
Because H has only real matrix elements in this basis, it is clear that that the secular 
det {z — H) = of such a system has only real coefficients, a fact closely linked to VT 
symmetry by the Bender-Mannheim theorem. 



III. SPECTRAL PROPERTIES OF PT-SYMMETRIC MODELS 

Spectral positivity plays a fundamental role in Hermitian systems. As we have seen, it 
ensures that the partition function Z is always positive, and it gives a representation of 
two-point functions as sums of decaying exponentials. On a practical level, it allows for the 
isolation of the lightest state in a given channel from the large-distance behavior of two- 
point functions. Because PT-symmetric models lead naturally to two different basis sets, 
completeness and the subsequent derivation of the Kallen-Lehmann representation for two- 
point functions is more complicated than in the Hermitian case. We give a self-contained 
derivation of both below, followed by a brief discussion of the implications of our results. Of 
particular interest is the connection of broken VT symmetry of the ground state (region III), 
with generalized Yang-Lee phase transitions. For simplicity, we will use H and /3 throughout 
this section, but we note that in classical statistical mechanics problems, the relevant objects 
are T and L, in which case /3 is not the inverse temperature of the system and correlation 
functions are between different spatial locations. As we will see for two-dimensional QCD 
in section |V], a transfer matrix approach can also be useful for quantum systems as well. 
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A. Completeness and Kallen- Lehman representation for "PT-symmetric systems 



We now prove a completeness relation for PT-symmetric models which is valid in all 
three regions. We consider the typical case where the Hamiltonian H is diagonalizable via a 
similarity transformation and has discrete, non- degenerate eigenvalues. Exceptional points, 
in the sense of degenerate real eigenvalues, occur at the boundary of region I as well as inside 
region II and III. There are well-known difficulties in the spectral resolution at these points; 
see [3j and [9j for specific examples as well as the general discussion in [lOj. The passage 
from region I to region III, where the ground state eigenvalue becomes degenerate, is a 
critical point analogous to those found in Hermitian models; the best-known example is the 
Yang-Lee singularity. In analogy with symmetry-breaking behavior in Hermitian models, 
it is likely necessary to define the behavior at such a point using either an infinitesimal 
symmetry-breaking perturbation or symmetry-breaking boundary conditions. At this stage 
of our understanding of PT-symmetric models in statistical mechanics, we cannot provide 
a general prescription, but must handle each model on a case-by-case basis. The case of 
non-degenerate eigenvalues, in contrast, is tractable for all models. By virtue of the secular 
equation for the eigenvalues, H and are isospectral: the existence of an eigenvalue- 
eigenvector pair for H 



(8) 



implies the existence of a corresponding pair for H . 



H"" |j> = E, I j> . 



(9) 



From the commutation relation [H, VT] = 0, we have 



HV \T]) = HVT |j) = VTH |j) = VTE, |j) = E*V \T]) 



(10) 



and its Hermitian conjugate 



{T3\VH^ = E,{T3\V. 



(11) 



Noting that VHV = H* we have 



(Til VH^ = (TjI V iVHVf = (TjI H^V 



(12) 



so that we find after multiplying by V on the right that 



{rj\H^ = E,{rj 



(13) 
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Now consider the matrix element 

Vj I {H^ - H"") I A;) = {E, - Ek) (rj\k) = (14) 
The set of eigenstates and | ^'^^ both complete so we must have 

(T]\k^^5,k (15) 

and also 

(tMj) oc (16) 
Thus we have two completeness relations 

and 

^Wf^'- '''' 

These completeness relations assume the Hamiltonian H acts in a space in which the inner 
products {n\Tn) are finite. The consequences of this assumption can vary from model 
to model. See [6j for a specific example where the contour along which a PT-symmetric 
differential equation is evaluated is deformed into the complex plane to ensure finite inner 
products. However, this technical point will not be of concern for the models we analyze 
in this paper. We will generally choose the normalization of states such that (Tj'lj) = 

{rj\j) = 1. 

With completeness relations at hand we can now calculate the partition function and 
spatial two-point functions in terms of eigenvalues and eigenstates. For a quantum theory, 
the partition function Z is given by 

Z = J2{n\ e"^^^^ \n) (19) 

n 

where {|n)} is an arbitrary orthonormal basis; note that the basis formed by the eigenstates 
of HpT cannot be used in this way to define Z unless HpT is Hermitian. When H is not 
Hermitian and PT symmetry is unbroken, it is a common practice to require VT\j) = \j) 
and to introduce an additional linear operator C such that CVT\j) • = 5jfc IH El- Using 
these conventions, one can write the partition function as a sum over eigenstates of the 
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form Z = (j| Ve~^^ \ j) where V = {CV)^ . The need for the operator V arises because 
the eigenstates \j) are not orthonormal with respect to the inner product {m\n). However, 
if H breaks PT-symmetry, not all of the eigenstates of H are eigenstates of PT, and an 
alternative expression for V of the form V = I'^j) would need to be used. We 

avoid these considerations in this paper by performing the trace in the partition function 
over an arbitrary orthonormal basis; such a basis may be constructed from the eigenstates 
\j) using the Gram-Schmidt algorithm. 
We insert the VT completeness relation 

^ - E ("I E ^ w - E ^"'"gif"^ - E (20) 

so that Z has the same form as in Hermitian theories, as assumed in section |Tlj We can 
write Z usefully as 

Z = J2 e-'^^" + + ^~^''') (21) 

r p 

where the sum over r is over all real energies and the sum over p is over pairs of complex 
energies. The oscillatory character of the second sum leads to negative contributions to 
the partition function, which is the sign problem. Note that in region I, there is no sign 
problem. The sign problem only arises in region II and III. Strictly speaking, the sign 
problem disappears in region II in the limit /3 — )■ oo. For transfer matrix problems, this 
implies that in the limit L — )■ oo, the sign problem becomes negligible in region II. Only in 
region III does the sign problem survive the infinite volume limit in the calculation of Z. 

B. Spectral theorems 

We now prove some general results for the typical case where H is symmetric but complex. 
When H is symmetric, we have \j) = |j) and \Tj) = and we immediately obtain the 
somewhat simpler completeness relation 

El^')^^^'^! (22) 

j 

In typical PT-symmetric quantum mechanics models, H is symmetric, and the completeness 
relation takes the form [U [3] 
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showing the exphcit normahzation factor. 

We begin by proving a result for certain matrix elements that appear repeatedly in these 
calculations. If Ej ^ E*, then the properly normahzed eigenstate of H with eigenvalue E* 
is related to VT \ j) by a phase factor: 

^r|j)=e-^|j*). (24) 

It immediately follows that 

VT \3*) = VTe-'^^Vr |j) = e'"^ |j) , (25) 

or 

gia,.* ^ gia, _ ^26) 

On the other hand, if Ej is real, we have simply 

^r|j) = e-^|j). (27) 

so in this case we may identify with We will prove results for matrix elements in the 
case where Ej is complex, and similar results will hold when Ej is real with replaced by 

li)- 

Consider the action of the bra 



{rf\vr 



(28) 



on an arbitrary ket 



XI \kc) + Xk* \K)) + ^Pn 



(29) 



pairs 



reals 



We have 



{rf\vr\^) = {Tf\vT 



X] (Xfe \kc) + Xk* \K)) + X] Pn 



jpairs 



reals 



{rf\vrm^{rr 



E (^^^'"' \K)+Xle'-' \kc)) + E p:e^"" \nr 



(30) 



(31) 



jpairs 



reals 
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and 



(32) 



(33) 



Now we apply the above result to a state of the form \ where VT(f)VT = 4>- 



{rk*\ci),\f) = {rk*\PT(p,PT\f) 



(34) 



{Tk* \(j>i\j*) ^ {Tk* iPme^'^lj) 



(35) 



(36) 



{Tk*\<i>^\f)^{rk\ci>,\j)* 



(37) 



This result has been proven for the case where both \j) and \k) are half of a conjugate pair 
of eigenvectors. It also holds when one or both of the eigenvectors has a real eigenvalue, in 
which case we identify \j*) with \j) . In the case where both Ej and Ek are real, the matrix 
element is real: 



With these results in hand, we now turn to one- and two-point functions in the case 
that H is symmetric and fields satisfy VT(f)VT — 0. We remind the reader again that 
in apphcation to classical statistical mechanics, the Hamiltonian is related to the transfer 
matrix by T = exp {—H), and /3 is identified with the length L of the system. A one-point 
function can be written as 



(rfc|</.i|j) = (rfc|0i|jr- 



(38) 



Zj 



jk 



jk 
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where we have used r to label the coordinate of x in the direction along which H acts. This 
expression for (0 (x)) is manifestly real as a consequence of our results for matrix elements. 

Similarly, suppose x and y are two points separated by a distance r in the direction in 
which H acts. Then the two point function for operator 0i (x) and 02 (z/)is defined by 

(01 (a;) 02(2/)) = ^Tr [0ie-'^^^-02e-(^-^)^^-] 

= i 5^e-'-^^e-(^-')^'= (Tfc |0i| j) (Tj |02| fc) (40) 

jk 

which is real using the same reasoning that was applied to the one-point function. This 
representation establishes the fundamental observable distinction between regions I, II and 
III. Each term in the spectral representation of the two-point function depends on r as 
exp [{Ek — Ej) r]. In regions I and II, the ground state is unique, and the terms with k = 
dominate for large /3. In region I, this leads to monotonic exponential decay. In region II, 
some of the excited states have complex energies, leading to modulated exponential decay 
in two-point functions. 

In region III, the ground state, defined as the state with the lowest value of Re (E), is no 
longer unique. The states |0) and |0*) will dominate in both Z and in two-point functions 
in the limit /3 — > oo, or in the limit L — t- oo for transfer matrices. We can take /3 sufficiently 
large that all states except Eq and Eq can be neglected, in which case Z can be approximated 
by 

Z ~ e-^^° + e-'^^o (41) 
and the approximate zeros of the partition function will occur at 

7„(i;„)=(3?!±l)l (42) 

where p is any integer. This is consistent with a general theory of partition function zeros 
that can be applied to models with PT-symmetric transfer matrices [llj. Under some 
technical conditions, the partition function in a periodic volume V = L'^ can be written as 

Z = J2 ^'^""^^ + O (e-^/^«e-^^^) (43) 

m 

where / = min^ Re [/m] and Lq is of the order of the largest correlation length of the system. 
The /m's have the interpretation of complex free energy densities, and are independent of L. 
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These phases are stable if Re (fm) — f or metastable otherwise. The zeros of the partition 
function are within O (e~^/^°) of the solutions of the equations. 

ReiU = Re{fn)^f 

Imifm) = Im{U) + {2p + l)^ 

for some m ^ n and p & Z . We can apply this directly to region III, using the representation 

reals pairs 

of the partition function where we identify /3 with L. We identify LEq and LEq as ^L'^fo 
and so that the partition function has a zero for values of the parameters such that 

Plm\h]=^^^^ (46) 

This tells us that the zeros of the partition function he on the boundary of region III, defined 
by Im [/o] = 0, in the limit V ^ oo. As the volume of the system is taken to infinity, the 
zeros of the partition function lie asymptotically on the boundary between phases. Note that 
this analysis depends on Lq remaining finite. At a 2nd-order transition, Lq goes to infinity 
and the approximation is invalid. Zeros of Z in region III can lead to potentially rapid 
oscillation of n-point functions , and correlation functions are ill-behaved in the vicinity of 
such points. 



IV. VT SYMMETRY IN CLASSICAL STATISTICAL MECHANICS 

In this section we consider several models of classical statistical mechanics that are VT- 
symmetric. In all of these, it is convenient to discuss the one-dimensional version of the 
model, which is analytically tractable. In addition to models where the classical Hamilto- 
nian, and hence the transfer matrix, is complex, there are also models where the classical 
Hamiltonian and transfer matrix are real, but the transfer matrix is not symmetric. Such 
models have a "hidden" VT symmetry. In those cases where the classical Hamiltonian is 
real, the matrix elements of T are positive. The Perron- Probenius theorem applies, and the 
eigenvalue of T of greatest magnitude will be real. Such models thus may lie in region I or 
II, but never in region III. 
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A. Z{N) Models 



Z{N) spin systems with complex magnetic fields arise naturally as simplified models of 
SU{N) gauge theories non-zero chemical potential [TB], with the case = 3 correspoinding 
to QCD at finite baryon density. These models are naturally PT-symmetric. In these 
models, there is a clear connection of VT symmetry with Z{N) Fourier transforms. This is 
not surprising: The Fourier transform of a real function f{x) on R obeys 

fikr = n-k) (46) 

or equivalently 

f{-kr = f{k) (47) 

which is precisely the statement of VT symmetry. In spin models with complex weights 
based on groups such as f/(l), Z{N) and SU{N), VT symmetry implies that the character 
expansion of a PT-symmetric model has real coefficients. This explains why the fiux-tube 
model [371 138] gives a purely real representation of the same physics of the Z{3) model to 
which it is dual. VT symmetry plays a similar role in the worldline approach to lattice field 
theories at non-zero chemical potential [14J . 

On each lattice site j of a Z{N) spin model there is a spin Wj, an element of the group 
Z{N) which may be parametrized as Wj = exp {2'Kinj/N) with rij G {0, 1, ...,N — 1} defined 
modulo so that and are identified. We take the operator V to be charge conjugation, 
acting as rij — )■ —rij, or equivalently Wj — )■ w*. The operator T is again complex conjugation. 
Although V and T have the same effect on the w/s, one is a linear operator and the other 
antilinear. We will show below that V is implemented as a unitary matrix in the transfer 
matrix formalism. The classical spin-model Hamiltonian H is defined by 

- = ^ ^ {wjwl + w*Wk) + J2 [^R (^i + ^j) + (^J - ^T)] (48) 

where /3 = 1/T, J, hji and hj are real and the sum over {jk) represents a sum over nearest- 
neighbor pairs. H is trivially PT-symmetric. This class of models has complex Boltzmann 
weights for > 3 when hj ^ 0. In the one-dimensional case, it is convenient to write H in 
the form 



(49) 
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where Hi and H2 are also real parameters. The partition function is given by the sum over 
all spin configurations 

Associated with the Hamiltonian is a transfer matrix such that Z — Tr T^' , where 



Tjk = exp 



^ {z^z'^* + z^*z^) + ^ {z^ + z'^) + +^ {z*^ + z*'') 



(51) 



and Ng is the spatial size of the lattice, and z — exp \2m/N\ is the generator of Z{N). The 
allowed values of j and k can be taken to run over either the set {0, 1, .., A?" — 1} or the set 
{1, 2, .., N}, and we generally identify the indices and N. 

We define the paritiy operator V by Vjk = Sj^^-k-, which satisfies — 1. Because 
z^~^ — z*^ , it is easy to see that 

VTV = T* (52) 

an equation also satisfied by the Hamiltonian when written in matrix form. This is the 
fundmental relation of VT symmetry, and can also be written as [VT, T] = 0. 
The discrete Fourier transform, defined by 

J'jk = ^z^', (53) 



is a symmetric and unitary operator satisfying J^T'^ — I. Furthermore, — T^"^ — V, so 
we have that the Fourier transform of T, T, obeys 

f* = {FTF+y 

— p+p*p 

= F'^PTPF 

= F+F^T{F+fF 
= FTF^ 
= f 

so we see that the Fourier transform of the transfer matrix is indeed real. 

The T'T-symmetric spin models have a complex order parameter coupled to a complex 
external field. We can prove that VT symmetry implies that the order parameter is real. 
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Let w he a spin operator in a VT -symmetric Z{N) spin system. The expectation value of 
w is given by 

(w) = ^ (54) 

where Z = Tr [T^"] and A^^ is the extent of the system in the direction in which the transfer 
matrix T acts. Using VT symmetry, we have 

Tr \PwP^T^'P] 
(w) = ^ - 

Tr [w*T*^^] 
^ Z 

= 

where we have used the fact that Z is real. The reality of (w) is analogous to the result 
that (x) is purely imaginary in PT-symmetric quantum mechanics. An alternative proof of 
the reality of (w) can be given: in the representation of the system induced by the discrete 
Fourier transform, the coefficients of the matrix representations of both w and T are all real. 
Equivalently, this can be also be seen easily from the character expansion of expressions like 

N-l 

exp [Hiz + H2Z*] = ^ ajz^ (55) 

j=0 

where all the coefficients aj are real if ifiand H2 are real. 

We illustrate the rich behavior possible in these models using the case of a Z{3) model 
in = 1 f3l]. If hj = , then the transfer matrix T is Hermitian. When hj 7^ 0, —/3H is no 
longer real and T is no longer Hermitian, but is VT symmetric. Figure [T] shows the phase 
diagram in the — hj plane for J = 0.2. There are four distinct regions. In region la, all 
three eigenvalues of the transfer matrix are real and positive. This region includes the line 
hi = 0, and has properties similar to those found in the Hermitian case. In region lb, all of 
the eigenvalues are real, but at least one of them is negative. In region II, the eigenvalue of T 
largest in magnitude is real, but the two other eigenvalues form a complex conjugate pair. In 
region HI, the two eigenvalues largest in magnitude form a complex conjugate pair, and the 
third, smaller, eigenvalue is real. In both region II and region HI, VT symmetry is broken, 
but in different ways. Borrowing the terminology from PT-symmetric quantum mechanics, 
we will describe the behavior in region HI as PT-symmetry breaking of the ground state, 
while region II is VT - symmetry breaking of an excited state. The behavior of the two-point 
function G {\j — k\) = (^w (j) (A;)) differs substantially in the three regions. In region I, 
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Figure 1: Phase diagram for the d = 1 "PT-symmetric Z{3) spin model in the — hj plane at 
J = 0.2. The interpretation of regions la, lb, 11 and 111 are given in the text. 

the two-point function falls off exponentially. We show typical behavior in region la in figure 
[2]for point A where {hji,hj) = (—0.45,0.5). Similar behavior occurs in region lb, as shown 
in the figure for point B where {h^, hj) = (—2.0, 1.5). Although the figure shows that the 
continuation of the two-point function away from integer values can be negative, note that 
the values at integer points are all non-negative. The two-point function at point C in region 
II where {hn, hj) = (0.25, 1.25) shows the damped oscillatory behavior associated with VT 
breaking in excited states. For the point D in region III, where {hn, hj) = (—0.5, 0.875), the 
VT breaking of the ground state leads to oscillatory behavior of the two-point function in 
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Figure 2: The two-point function as a function of lattice spacing for the parameters corresponding 
to points A, B, C and D in figure [T] 

the limit of large distance. Note that region III only occurs when hn is negative. For < 
and hi = 0, the spin configurations with lowest energy have a two-fold degeneracy. With 
hi = 0, the ground state of the transfer matrix is unique. For the case h^ < 0, hi = 0, and 
J large, the splitting of the two lowest eigenvalues of the transfer matrix in = 1 is small. 
For sufficiently strong hi, the real parts of the two lowest eigenvalues of T merge, and VT 
symmetry breaking of the ground state occurs. 



B. The chiral Potts model 



Out first example of a system with a hidden VT symmetry is the chiral Potts model 
IB] . It is a variant of Z{N) spin models, and we use the same notation as above. 
Consider a d = 1 Z{N) spin model with Hamiltonian "H of the form 

- /3'H = ^ ^ [wjuw*_^_i + w*u*Wj+i) (56) 
j 

where u = exp {27TiA/N) with A G [0,1] . The classical Hamiltonian "H is real, and the 
model is invariant under the action of T, regarded as complex conjugation. However, the 
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transfer matrix is not Hermitian for general A. The transfer matrix is given by 

J 



Tjk = exp 



(57) 



where as before z = exp (27Ti/N). The reahty of "H imphes that T*^ = Tjk, and thus possesses 
a generahzed VT symmetry, with V simply taken to be the identity. Although T is real, 
it is not symmetric in general, because 7^ T. Because the transfer matrix has only 
real positive entries, by the Perron-Frobenius theorem the correlation functions of the chiral 
Potts model can exhibit region I or region II behavior, but never region III. 

(58) 

C. The ANNNI model 



The Anisotropic Next-Nearest-Neighbor Ising (ANNNI) model [39j is a prototypical ex- 
ample of a system that appears as if it "should" have a Hermitian transfer matrix, but does 
not. Instead, the model has a generalized VT symmetry which underlies the model's un- 
usual phase structure [31] • The one-dimensional model is exactly solvable, and has a reduced 
Hamiltonian 

/3H = -Ki ^ SjSj+i - ^ SjSj+2 (59) 
j j 

where Ki and K2 are real couplings and the Ising spins take on the values ±1. The Hamil- 
tonian is real, and stochastic simulations of the model may be carried out with ease. One 
approach to solving the model is to construct a 4 x 4 transfer matrix between nearest neighbor 
pairs 



/ 2K1+2K2 



V 



-2K2 



^2K2-2Ki g-2i^2 
g-2X2 ^2K2-2Ki 



^2Ki+2K2 



(60) 



e ^ e ^ e ' e 
The partition function for N spins, with even and periodic boundary conditions, is given 



hj Z = Tr 



N/2 



. The matrix T4 is real but not symmetric. It commutes with a generalized 
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V 



(61) 



parity operator V of the form 

/O l\ 
10 
10 

\ 1 0/ 

which implements the symmetry of the model under s — )> — s; T is complex conjugation, and 
acts trivially on T4. The combination of the reality of T4 together with [T4, P] — Q establishes 
that T4 is PT-symmetric. There is another approach to solving the one-dimensional model 
which better displays its VT symmetry. We introduce a set of bond variables Oj into 
the partition function Z which we force to be equal to SjSj+i via the Z2 delta function 
(1 + ajSjSj^i) /2. We can then write = CTjCrj+i in the Hamiltonian. It appears that 

the new Hamiltonian is simply 

j 3 

and the model reduces to a standard Ising model in an external field. This is somewhat 
misleading, because there remains a determinantal factor associated with the change of 
variables 



EEn 



l + ajSjSj+i 



(63) 



Carrying out the sum over the s variables with periodic boundary conditions, we find 



i+n<^. 



(64) 



which tells us that there is a global constraint on the partition function: only configurations 
with Y\j — 1 contribute. Let T2 be the 2x2 transfer matrix of the one-dimensional Ising 
model in an external field 



To 



(65) 



We define another matrix, T2, as 



1 
i 




(66) 



such that the transfer matrix of the model is the 4x4 matrix 7^ = T^j ©T^. The square of the 



eigenvalues of are the eigenvalues of Tlj, as they must be, and Z — Tr 



-^4 



The transfer 
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matrix T4 is invariant under PT, with the parity operator V given by 1 ©(13, and T given by 
complex conjugation. This construction leads directly to the same eigenvalues found in |39] . 
The eigenvalues of T2 are of course always real, while the eigenvalues of T2 are either real or 
form a conjugate pair. For coshi^i > e~^^^, the eigenvalues of T2 are real, and the spin-spin 
two-point function decays exponentially. The system is in region I. For coshi^'i < e~^^^, 
the eigenvalues of T2 are complex, and the spin-spin two-point function shows a periodic 
modulation of its exponential decay. This is region II, and the line coshi^'i = e"^^^ defines 
the disorder line separating the two regions. The eigenvalues of T2 are always smaller in 
absolute value than the eigenvalues of T2 in this model, so region III does not occur in the 
d = 1 ANNNI model. 

V. Vr SYMMETRY IN QUANTUM STATISTICAL MECHANICS MODELS 

All quantum many-body problems involving a non-zero chemical potential may be de- 
scribed in terms of a non-Hermitian Hamiltonian with generalized VT symmetry j3l]. At 
first glance, this is surprising, but it is a simple consequence of the use of Wick rotation and 
the Euclidean formalism for equilibrium statistical mechanics. This PT-symmetric descrip- 
tion is closely related to the sign problem. We will explain in detail how the sign problem 
arised in QCD with heavy quarks at non-zero chemical potential. The two-dimensional 
case will be solved numerically as an application of VT symmetry to this class of problems 
[331 EH- 

A. VT symmetry at finite density 

We start from a theory with a Hermitian Hamiltonian H and a conserved global quantum 
number N, obtained from a conserved current j'^, that commutes with H. We assume that 
H is Hermitian and invariant under the combined action of time reversal T and a charge 
conjugation C that reverses the sign of j". We take the number of spatial dimensions to 
he d — 1, and the spatial volume to be L'^^^. The grand canonical partition function at 
temperature T = and chemical potential fi is given hj Z = Tr [exp {—^H + fSfiN)]. If 
Z is written as a Euclidean path integral, the time component of the current j° will Wick 
rotate to ij'^, while the chemical potential fi does not change. This leads directly to a non- 
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positive weight in the path integral, and is the origin of the sign problem in finite density 
calculations. The Euclidean space Lagrangian density may be written as £ — i/ij'^ where 
jC is the Euclidean Lagrangian for = 0; £ — ifij'^ is complex. The nature of the problem 
is changed by changing the direction of Euclidean time, so that we are now considering a 
problem at zero temperature with one compact spatial dimension of circumference f3. Upon 
returning to Minkowski space, p does not rotate. We pick, say, the u — 1 direction to be 
the new time direction and the new inverse temperature is L. When /i — 0, the original 
Hamiltonian is obtained, but for // 7^ the partition function is now given by 

Z^Tr [e-^^p] (67) 

where 

Hp^H-i/i J d'^-^xj''. (68) 

The new Hamiltonian Hf^ is non-Hermitian, but possesses a generalized VT symmetry, where 
the role of V is played by the charge conjugation operator C that changes the sign of j° and 
N. Under the combined action of CT, j'^ — >■ —j'^ and i — > —i, leaving the Hamiltonian HpT 
invariant. If we introduce the operator Hl — H — /iN, we have the relation 

Z = Tr [e-"''^] = Tr [e"^^"] (69) 

induced by the space-time transformation that exchanges directions 1 and d. Note that Z 
is obtained from Hj^ by a sum over all eigenstates, but is dominated by the ground state of 
Hp in the limit of large L. 

B. d = 2 gauge theories 

Within the Euclidean space formalism, a non-zero temperature T is obtained by making 
the bosonic fields periodic in Euclidean time, with period /3 = 1/T. On the other hand, a 
non-zero chemical potential must be implemented in a way that makes the weight function 
used in the Feynman path integral complex, as we have seen above. We will show below 
exactly how QCD with quarks at finite density may be interpreted as a theory with VT 
symmetry. 

The Polyakov loop plays a crucial role. Defined as a path-ordered exponential of the 



25 



gauge field, in 3 + 1 dimensions the Polyakov loop operator P is given by 



P(x) = Pexp i dtAi{x,t) , 



(70) 







and represents the insertion of a static quark into a thermal system of gauge fields at a 
temperature T = Figure [s] shows the Polyakov loop in this geometry. Because of 

the periodic boundary conditions in the Euclidean time direction, the Polyakov loop is a 
closed loop, and its trace is gauge invariant. Also known as the Wilson line, the Polyakov 
loop represents the insertion of a static quark at a spatial point a; in a gauge theory at 
finite temperature. In particular, the thermal average of the trace of P in an irreducible 
representation R of the gauge group is associated with the additional free energy required 
to insert a static quark in the fundamental representation via 



Pure SU{N) gauge theories have a global Z{N) symmetry P — )■ zP where z = is the 
generator of Z{N), the center of SU{N). This symmetry, if unbroken, guarantees that for 
the fundamental representation F, {TrpP (x)) = 0. This is interpreted as Fp being infinite, 
and an infinite free energy is required to insert a heavy quark into the system. On the other 
hand, if the Z{N) symmetry is spontaneously broken, the free energy required is finite. 
Thus confinement in pure gauge theories is associated with unbroken center symmetry, and 
broken symmetry with a deconfined phase. The Polyakov loop is the order parameter for 
the deconfinement transition in pure gauge theories (TrpP) = in the confined phase and 
(TrpP) 7^ in the deconfined phase. The addition of dynamical quarks in the fundamental 
representation explicitly breaks this Z{N) symmetry. Nevertheless, the Polyakov loop re- 
mains important in describing the behavior of the system, as we will see in our treatment 
of the sign problem. 

In pure gauge theories, the Wilson loop operator is used to measure the string tension 
between quarks in the confined phase where vanishes for representations transforming 
non-trivially under Z{N). At non-zero temperature, a timelike string tension cr[*'' between 
k quarks and k antiquarks can be measured from the behavior of the correlation function 




(71) 



( 



TrpP'' {x)TrF (P+ (y))') 



~ exp 




x-y 



(72) 
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Figure 3: The Polyakov loop in Euclidean space-time. 

at sufficiently large distances. A confining phase is defined by two properties: the expectation 
value (Tr/jP) is zero for all representations R transforming non-trivially under Z{N), and 
the string tensions o"^*'' must be non-zero for A; = 1 to — 1. 

Perturbation theory can be used to calculate the one-loop free energy density fq of quarks 
in (i + 1 dimensions in the fundamental representation with spin degeneracy s moving in a 
Polyakov loop background at non-zero temperature T = and chemical potential /i 

= -sT [ -^^Trn [in (l + Pe^^-^^'^) + In (l + p+e-^f'-^'^^)] (73) 

where cok = y/k"^ + M"^ is the energy of the particle as a function of k and M is the mass 
of the particle |2T| W2\ . The expression for a bosonic field is similar. The logarithm can be 
expanded to give 

This expresssion has a simple interpretation as a sum of paths winding around the timelike 
direction. With standard boundary conditions, which are periodic for bosons and antiperi- 
odic for fermions, this one-loop free energy always favors the deconfined phase. 

The effects of heavy quarks in the fundamental representation, with /3M ^ 1, on the 
gauge theory can be obtained approximately from the n = 1 term in the free energy 

/, ^ -sT [ -^Trp [Pe^'^-^^' + p+e"^^"/'-^] (75) 
J (27r) 
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Figure 4: The Polyakov loop in (1 + l)-dimensional space-time. 

because term with higher n are suppressed by a factor e""^*^. In this approximation, bosons 
and fermions have the same effect at leading order. After integrating over /c, the free energy 
fq can be written as fg —hp \e^^^TrpP + e'^^TrpP^']^ . The one- loop free energy density 
is the one-loop effective potential at finite temperature. Thus the free energy for the heavy 
quarks can be added to the usual gauge action to give an effective action which involves 
only the gauge fields. The effective action is given by 



I d^+^x - hp (e^^'TrpiP) + e-^^'Trp{P+)) 



(76) 



and the structure and symmetries of the theory are obviously the same in any number of 
spatial dimensions. Because TrpP is complex for > 3, the effective action for the gauge 
fields is complex. This is a form of the so-called sign problem for gauge theories at finite 
density: the Euclidean path integral involve complex weights. This problem is a fundamental 
barrier to lattice simulations of QCD at finite density. 

C. Heavy quarks at /x 7^ in two dimensions and T'Tsymmetry 

In one space and one time dimension, the field theory arising from the effective action 
Seff can be reduced to a PT-symmetric Hamiltonian acting on class functions of the gauge 
group. The effective action, including the effects of heavy quarks, is 

Seff = jd'^ ^2 (^;.)' - {e^'Trp{P) + e-^^Trp{P+)) (77) 

where the gauge field now has two components. Figure |4] shows the Polyakov loop in a 
1 -|- 1-dimensional geometry. It is convenient to work in a gauge where Ai = 0; this is turn 
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Figure 5: The Polyakov loop in a (1 + l)-dimensional transfer matrix geometry, 
implies that A2 depends only on Xi. After integration over X2, we are left with a Lagrangian 

L = ^, (^) ' - [e^'TMP) + e-^^rr^(P+)] (78) 

which we regard as the Lagrangian for a system evolving as a function of a time coordinate 
Xi. This represents a change from a Euclidean time point of view to a transfer matrix 
geometry, as shown in figure [5| In this geometry, the Polyakov loop represents the insertion 
of an electric flux line in a box with periodic boundary conditions, and the free energy 
density is obtained from the lowest-lying eigenvalue of the transfer matrix. 

The physical states of the system are gauge-invariant, meaning that they are class func- 
tions of P: \E' [P] = \1' [gPg^]. The group characters form an orthonormal basis on the 
physical Hilbert space: [P] = 'Y^^a^iTrji (P). The Hamiltonian obtained from L, acts 
on the physical states as 

H = - hpP [e^^TrpiP) + e-^^Tr^(P+)] (79) 

where C2 is the quadratic Casimir operator for the gauge group, the Laplace-Beltrami op- 
erator on the group manifold. We have thus reduced the problem of heavy quarks at finite 
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density in two dimensions to a problem of quantum mechanics on the gauge group. Unfor- 
turnately, the Hamiltonian H is not Hermitian when /i 7^ 0, and thus cannot be rehed upon 
to have real eigenvalues. This is a direct manifestation of the sign problem. 

Although the Hamiltonian H is not Hermitian when /i 7^ 0, it is "PT-symmetric under 
the transformations 

V:x2^ -X2 A2 -A2 (80) 



r-.i^-i (81) 

which should be regarded as parity and time-reflection in the transfer matrix geometry. 
Together these lead to 

VT:P^P (82) 

which leaves the Hamiltonian invariant. If this PT-symmetry is unbroken, the eigenvalues 
of the Hamiltonian will be real, and there is no sign problem. The VT symmetry remains 
even in the high-density limit where the quark mass M and chemical potential /i are taken 
to infinity in such a way that antiparticles are suppressed and P+ does not appear in H. 

The simplest non-trivial gauge group is SU{'i), because the cases of U{1) and SU{2) are 
atypical. For the gauge group f/(l), the Hamiltonian H may be written as 

H = -'-Y^2-h^^ (e^^+^^ + e-^^-^') (83) 
but a simple change of variable 6^6 + if3fj, eliminates fi: 

H=-"~2-hFn^^''+^-'') m 

This is very similar to the case of the two-dimensional PT-symmetric sine-Gordon model 
considered in [8]. In the case of SU(2), all the irreducible representations are real, and the 
Hamiltonian is Hermitian: 

Hsu(2) = ^C2 - 2hF cosh (/3/i) x,=i/2(P). (85) 

In general there is no sign problem in SU (2) gauge theories at finite density holds in general, 
and this feature has been exploited in lattice simulations with 7^ |23l 126] . 

Thus = 3 is the first non-trivial case for SU{N) gauge groups. We have calculated the 
lowest eigenvalues of H using finite dimensional approximants. It is convenient to work in 
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Figure 6: The real part of the SU{3) Hamiltonian Hj^ as a function of Pfi. The upper graph is 
for periodic boundary conditions for the heavy quarks, while the lower graph is for antiperiodic 
boundary conditions. The energy has been scaled such that is set equal to 1. 

the group character basis. The Casimir operator C2 is diagonal in this basis, and characters 
act as raising and lowering operators. For example, in the 4x4 subspace spanned by the 1, 
3 , 3, and 8 representations of SU{3), the Hamiltonian takes the form 

/ e-'^'^hp/S e'^'^hp/S \ 



4 



(86) 



If hp is set to zero, we see that the eigenvalues are proportional to Casimir invariants 0, 
4/3 , 4/3, and 3 for the 1, 3, 3, and 8 representations of SU{?>). We have therefore removed 
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an overall factor of so the overall strength of the potential term is controlled by the 

dimensionless parameter Ihpjg^. The resulting dimensionless energy eigenvalues are thus 
normalized to give the quadratic Casimir operator when 2hp/ g"^ = 0. The lowest eigenvalues 
have been calculated numerically using a basis of dimension nine or larger, with the stability 
of the lowest eigenvalues checked by changing the basis size. 

The parameter hp is positive for fermions with antiperiodic boundary conditions in the 
timelike direction, which are required for spectral positivity. However, it is also of interest 
to consider the case of periodic boundary conditions for the heavy quarks, corresponding to 
hp < PT| 1221 SI]- In figure |6| we show the real part of the eigenvalues of Hp, measured 
in units where is set to 1. The overall strength of the potential term is set by the 

dimensionless parameter 2hp/g'^. In the upper graph, 2hp/g'^ = —0.5, corresponding to 
periodic boundary conditions for the heavy quarks. The lower graph shows the real part of 
the energy eigenvalues for 2hp/ g"^ = 0.5. In both cases, we see the real parts of pairs of energy 
eigenvalues coalescing as 13 fi is increased. At the point where the real parts become identical, 
these energy eigenvalues acquire an imaginary part, indicative of broken VT symmetry. In 
the case of periodic boundary conditions, we see that the ground state shows PT-symmetry 
breaking before any of the higher states; thus for large this places the system in region 
III. Note that for > 3, the heavy quark finite density problem of SU (N) gauge theory is 
in the universality class of the Lee- Yang problem for Z{N) spin systems [16]. In the physical 
case of antiperiodic boundary conditions, "PT-symmetry breaking appears to occur only in 
excited states. In the case where all eigenvalues are real, which appears to hold for small 
PTsymmetry is unbroken and the system is in region I. This in turn implies that for small 

the sign problem may be solved in principle by a similarity transform to a Hermitian 
Hamiltonian. For large Pfi, the VT symmetry is broken in some of the excited states, which 
will lead to the region II behavior of sinusoidal decay of spatial correlation functions at high 
density. 



VI. CONCLUSIONS 



For Hermitian systems, there is a well-developed understanding of critical behavior and 
phase structure, connecting a wide range of systems from simple classical spin systems to 
exotic quantum field theories. For PT-symmetric models, we are in a sense starting over 
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again with a richer, larger class of systems. There are indications that VT symmetry is 
crucial to the understanding of the sign problem. More generally, the characterization of 
the phase structure and universality classes of T'T-symmetric systems is a logical extension 
of the successful effort to characterize critical phenomena in Hermitian systems. 

Acknowledgments 

Michael Ogilvie wishes to acknowledge the support of his research by the US Department 
of Energy. 



[1] Mark G. Alford, Andreas Schmitt, Krishna Rajagopal, and Thomas Schafer. Color supercon- 
ductivity in dense quark matter. Rev.Mod.Phys., 80:1455-1515, 2008. 

[2] A.A. Andrianov. Dual oscillators and Quantum Pendulums: spectrum and correlators. 
Phys.Rev., D76:025003, 2007. 

[3] A.A. Andrianov, F. Cannata, and A.V. Sokolov. Spectral singularities for Non-Hermitian 
one-dimensional Hamiltonians: Puzzles with resolution of identity. J.Math.Phys., 51:052104, 
2010. 

[4] Carl M. Bender. Introduction to PT-Symmetric Quantum Theory. Contemp.Phys., 46:277- 
292, 2005. 

[5] Carl M. Bender. Making sense of non-Hermitian Hamiltonians. Rept.Prog.Phys., 70:947, 2007. 
[6] Carl M. Bender and Stefan Boettcher. Real Spectra in Non-Hermitian Hamiltonians Having 

PT Symmetry. Phys. Rev. Lett, 80:5243-5246, 1998. 
[7] Carl M. Bender, Dorje C. Brody, Jun-Hua Chen, Hugh F. Jones, Kimball A. Milton, et al. 

Equivalence of a Complex PT-Symmetric Quartic Hamiltonian and a Hermitian Quartic Hamil- 

tonian with an Anomaly. Phys. Rev., D74:025016, 2006. 
[8] Carl M. Bender, H.F. Jones, and R.J. Rivers. Dual PT-symmetric quantum field theories. 

Phys. Lett, B625:333-340, 2005. 
[9] Carl M. Bender and Philip D. Mannheim. Exactly solvable T'T-symmetric hamiltonian having 

no hermitian counterpart. Phys. Rev. D, 78:025022, Jul 2008. 



33 



[10] Carl M. Bender and Philip D. Mannheim. PT symmetry and necessary and sufficient conditions 
for the reality of energy eigenvalues. Phys.Lett., A374: 1616-1620, 2010. 

[11] M. Biskup, C. Borgs, J. T. Chayes, L. J. Kleinwaks, and R. Kotecky. General theory of lee-yang 
zeros in models with first-order phase transitions. Phys. Rev. Lett., 84:4794-4797, May 2000. 

[12] John L. Cardy. Conformal Invariance and the Yang-Lee Edge Singularity in Two-Dimensions. 
Phys.Rev.Lett, 54:1354-1356, 1985. 

[13] John L. Cardy and G. Mussardo. S Matrix of the Yang-Lee Edge Singularity in Two- 
Dimensions. Phys.Lett, B225:275, 1989. 

[14] Shailesh Chandrasekharan. A New computational approach to lattice quantum field theories. 
PoS, LATTICE2008:003, 2008. 

[15] J. Cox, C. Gattringer, K. Holland, B. Scarlet, and U.J. Wiese. Meron cluster solution of 
fermion and other sign problems. Nucl.Phys.Proc.Suppl., 83:777-791, 2000. 

[16] Thomas A. DeGrand and Carleton E. DeTar. Phase Structure of QCD at High Temperature 
with Massive Quarks and Finite Quark Density: a Z(3) Paradigm. Nucl.Phys., B225:590, 1983. 

[17] Patrick Dorey, Clare Dunning, and Roberto Tateo. Spectral equivalences, Bethe Ansatz equa- 
tions, and reality properties in PT-symmetric quantum mechanics. J. Phys. A, A34:5679-5704, 
2001. 

[18] Patrick Dorey, Clare Dunning, and Roberto Tateo. The ODE/IM Correspondence. J. Phys. A, 
A40:R205, 2007. 

[19] Patrick Dorey, Clare Dunning, and Roberto Tateo. Prom PT-symmetric quantum mechanics 

to conformal field theory. Pramana, 73:217-239, 2009. 
[20] M.E. Fisher. Yang-Lee Edge Singularity and phi**3 Field Theory. Phys. Rev. Lett, 40:1610- 

1613, 1978. 

[21] David J. Gross, Robert D. Pisarski, and Laurence G. Yaffe. QCD and Instantons at Finite 

Temperature. Rev. Mod. Phys., 53:43, 1981. 
[22] Fritz. Haake. Quantum signatures of chaos. Springer- Ver lag, Berlin; New York, 2010. 
[23] Simon Hands, John B. Kogut, Maria-Paola Lombardo, and Susan E. Morrison. Symmetries and 

spectrum of SU(2) lattice gauge theory at finite chemical potential. Nucl.Phys., B558:327-346, 

1999. 

[24] S. Howes, L.p. Kadanoff, and M. Den Nijs. Quantum Model For Commensurate - Incommen- 
surate Transitions. Nucl.Phys., B215:169-208, 1983. 

34 



[25] H.F. Jones, J. Mateo, and R.J. Rivers. On the Path-Integral Derivation of the Anomaly for 
the Hermitian Equivalent of the Complex PT-Symmetric Quartic Hamiltonian. Phys.Rev., 
D74: 125022, 2006. 

[26] J.B. Kogut, D.K. Sinclair, S.J. Hands, and S.E. Morrison. Two color QCD at nonzero quark 

number density. Phys.Rev., D64:094505, 2001. 
[27] T.D. Lee and Chen-Ning Yang. Statistical theory of equations of state and phase transitions. 

2. Lattice gas and Ising model. Phys.Rev., 87:410-419, 1952. 
[28] E.Y. Loh, J.E. Gubernatis, R.T. Scalettar, S.R. White, D.J. Scalapino, et al. Sign problem in 

the numerical simulation of many-clcctron systems. Phys.Rev., B41:9301-9307, 1990. 
[29] M. P. Lombardo. QCD at Non-Zero Density : Lattice Results. J. Phys., G35:104019, 2008. 
[30] AH Mostafazadeh. Exact PT symmetry is equivalent to Hermiticity. J. Phys. A, A36:7081-7092, 

2003. 

[31] Joyce C. Myers and Michael C. Ogilvie. New phases of SU(3) and SU(4) at finite temperature. 
Phys.Rev., D77:125030, 2008. 

[32] Joyce C. Myers and Michael C. Ogilvie. Phase diagrams of SU(N) gauge theories with fermions 
in various representations. JHEP, 0907:095, 2009. 

[33] Michael Ogilvie and Peter Meisinger. VT Symmetry and QCD: Finite Temperature and 
Density. SIGMA, 5:047, 14 p., 2009. 

[34] Michael Ogilvie, Peter Meisinger, and Timothy Wiser. "VT" symmetry in statistical mechan- 
ics and the sign problem. International Journal of Theoretical Physics, 50:1042-1051, 2011. 
10.1007/sl0773-010-0626-5. 

[35] Michael C. Ogilvie and Peter N. Meisinger. PT symmetry and large-N models. J. Phys. A, 
A41:244021, 2008. 

[36] S. Ostlund. Incommensurate and commensurate phases in asymmetric clock models. Phys.Rev., 
B24:398-405, 1981. 

[37] Apoorva Patel. A Flux Tube Model of the Finite Temperature Deconfining Transition in QCD. 

Nucl.Phys., B243:411, 1984. 
[38] Apoorva Patel. More on the Flux Tube Model of the Deconfining Transition. Phys. Lett., 

B139:394, 1984. 

[39] W. Selke. The ANNNI model: Theoretical analysis and experimental application. Phys. Rep., 
170:213-264, 1988. 



[40] M.A. Stephanov. QCD phase diagram: An Overview. PoS, LAT2006:024, 2006. 

[41] Mithat Unsal and Laurence G. Yaffe. Center-stabilized Yang-Mills theory: Confinement and 
large N volume independence. Phys.Rev., D78:065035, 2008. 

[42] Nathan Weiss. The Effective Potential for the Order Parameter of Gauge Theories at Finite 
Temperature. Phys.Rev., D24:475, 1981. 

[43] Chen-Ning Yang and T.D. Lee. Statistical theory of equations of state and phase transitions. 
1. Theory of condensation. Phys.Rev., 87:404-409, 1952. 

[44] V.P. Yurov and A.B. Zamolodchikov. Truncated Conformal Space Approach to Scaling Lee- 
Yang Model. Int.J.Mod.Phys., A5:3221 3246, 1990. 

[45] A.B. Zamolodchikov. Thermodynamics Bethe Ansatz in Relativistic Models: Scaling Three 
State Potts and Lee- Yang Models. Nucl.Phys., B342:695-720, 1990. 



36 



